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Motivation 


v 


•  Particle-solidification  front 
interactions  play  a  crucial  role  in 
predicting  microstructures  of  metal- 
matrix  composites  (MMCs) 


MMCs:  Micron-sized  particles  embedded  in  a 
metal  matrix  -  acts  as  reinforcement 


Industies:  Automotive,  aerospace 
Properties:  lightweight,  good  wear  resistance 


Processing  of  MMCs 


Directionally  solidified  AI-17Cu/SiC.  21x 


As-cast  AI-2Mg/20  wt  pet  graphite.  31  x 


Want:  Uniform  distribution  of  particles 
Need:  Understand  particle-front  interactions 


What  happens  when  a  solidification  front  interacts 
^ with  a  particle? 


•  As  solidifying  front  approaches  a  particle  it  will  either  be 

-  Engulfed 

-  Pushed 

-  Pushed,  then  engulfed 


•  Critical  velocity,  VCR 

-  Solidification  velocity,  Vs  <  VCR  ->  particle  pushing 

-  Vs  >  VCR  ->  particle  engulfment 

♦  Challenge:  A  multiscale  model  that  couples  the  dynamics  in  the  gap  with  the 


overall  thermal/fluid  transport 


Physics  of  particle-solidification  front 
interactions 


Dynamics  of  the  particle 


•  Force  balance  on  the  particle 


dV 


™  - - -  =  Fi  ~  fd 


dt 


•  F,  is  the  repulsive  interaction 
force 

•  Fd  is  an  attractive  drag  force 

•  Balance  determines  whether 
the  particle  will  be 
engulfed/pushed 


Fluid  mechanics  and  thermal  transport 


•  In  addition,  thermal  conditions  affect  the  critical  velocity.  Large  scale  effects 
(thermal  conductivity  differences  between  the  particle  and  the  solidification  front), 
and  small  scale  effects  (premelting  film). 


The  thermal  conductivity  issue 


♦  kp/k,  has  a  significant  impact  on  whether  a  particle  is  pushed/engulfed 


♦  In  fact,  experimental  observations  have  shown  that: 

♦  Engulfment  1 00%  of  the  time  for  kp/k,  >  1 .0 

♦  Pushing  for  kp/k,  <  1.0 


♦Why  does  the  kp/k,  have  such  a  significant  impact? 


How  does  the  kp/k,  influence  the  interface  shape? 


kp/k,<  1 .0 


kp/k,  >1.0 


Summary:  The  thermal  conductivity  ratio,  kp/k|,  has  a  significant  impact 

on  whether  a  particle  is  pushed/engulfed  ->  this  is  a  micro-scale 
phenemenon  (scale  of  the  particle) 


What  about  nano-scale  phenomena? 


Premelted  films 
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Film  thickness:  d  =  (-2cr2 — — — )1/3r  1/3 

Pl*H/ 


Atomic  radius:  cr 


Liquid  density:  Pi 


Latent  heat  per  molecule:  AH 
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A r  =  rsp-(r,p+r,s ) 


Premelting  film  theory  effect  on  interface 
^temperature  of  a soljdification  front 


T,= 
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f-T 
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Tj  is  the  interface  temperature 
Tm  is  the  melting  temperature 


v=  3  for  van  der  Waals 

A  is  an  interaction  length  scale  (nanometer  scale) 


Disjoining  pressure  theory 


•  Disjoining  pressure  - ► 

-  Nano-scale  force 

-  Molecular  in  nature 

A  -  Hamaker  constant  ~A y 

Ar  =  rsp-(rip  +  ris) - 


A 

6nd 3 


Ay  >  0  ->  repulsive 
Ay  <  0  ->  attractive 


Van  der  Waals? 


Previous  work  on  thermal  effects  of  premelted  layers  on  particle- 

solidification  front  dynamics 

(Garvin  and  Udaykumar,  Journal  of  Crystal  Growth  (2005)) 


Premelting  film  theory  effect  on  interface 
temperature  of  a  solidification  front 


Vpull  =  1000  micron5/5ec  -  includes  premelting  effects 


Vpull  =  1000  microns/sec  -  No  premelting  effects 
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The  premelting  effect  on  VCR 


2  approaches  when  modeling 

1)  Not  to  include  premelting  in  the  interaction 

Requires  setting  a  gap  thickness  cut-off  value  for  engulfment 
(ad  hoc) 

2)  Including  premelting  (Rempel  and  Worster) 

-  Linear  stability  analysis,  steady-state  and  k  /k,  =  1 .0 

-  No  need  for  ad-hoc  cut-off  for  gap  thickness 


Summary  of  physics 


kp/k|  ratio  greatly  affects  the  critical  velocity  (micro-scale  effect) 
-  Alters  the  interface  shape 

•  kp/k|  <  1 .0,  bump  forms  that  promotes  pushing 

•  kp/k|  >  1 .0,  trough  forms  that  promotes  engulfing  (higher 
drag) 


Premelting  film 

-  Liquid  layer  wants  to  exist  between  the  front  and  the  particle 

•  Also  alters  interface  shape 

•  If  no  premelting  is  assumed,  gap  thickness  cut-off  value  is 
needed 

The  kp/k,  effects  act  on  a  larger  (microscale)  whereas  the 
premelting/disjoining  pressure  effects  act  on  a  “nano”-scale 


PROBLEM:  How  to  model  particle-solidification 

front  interactions? 


In  the  past ... 


-  Theoretical  research  relied  on... 

1)  All-simplified  models  for  the  forces  acting  on  the  particle 

2)  All-qualitative  explanations/models  based  off  of  thermal 
property  considerations 

3)  Most  on  steady-state  analysis 

4)  Most-  Do  not  assume  a  premelting  layer  in  the  gap 

5)  All-assumed  stable  (non-dendritic  solidifications  front) 

ALL  MODELS  ATTEMPT  TO  FIND  CRITICAL  VELOCITY 


OUR  APPROACH... 

Lubrication  Model  for  the  flow  in  the  gap  between  the 
front  and  particle  and  NS  for  outer  flow  -  match 
solutions  (matched  numerics) 


Advantages 


1.  Capture  dynamics  (i.e.  not  a  steady-state  solution) 

2.  Premelting  effects  will  be  included 

3.  kp/k,  effects  will  be  captured  from  the  governing  equation 

4.  Do  not  have  to  rely  on  simple  force  models  as  the  lubrication 
equation  will  contain  a  repulsive  body  force 

5.  Can  be  generalized  for  a  variety  of  problems  involving  solid 
body  interactions  (including  dendrite+particle) 


Modeling  of  particle-front  interactions 


Main  Ingredients: 

1 .  Solidification  front  and  the  particle  are  tracked  using  a  sharp-interface  level- 
set  method  with  transport  variables  solved  via  NS  equations  “outer  region” 

2.  Solution  in  the  gap  in  treated  as  a  lubrication  layer,  so  a  different  set  of 
equations  are  needed  “inner  region” 

3.  The  solution  of  the  “inner”  and  “outer”  regions  are  coupled  through  BCs. 

“matching  region” 


liquid  particle 


solid 


Outer  transport  equations 


•  Navier-Stokes  Equations  (with  mass  and  energy  conservation): 


Continuity  equation: 


du  dv 
—  +  —  =  0 
dx  dy 


x-Momentum  equation: 


y-Momentum  equation: 
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Energy  equation: 
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Interface  conditions: 


pHk^}=K 
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Level-Set  Method 


Level-set  function  (§)  -  Signed  normal  distance  from  the  interface 


Bulk  Points 


Interfacial  Points 
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(O-levelset) 
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Curvature 


Using  a  Lubrication  Model  for  “inner”  region 


•  Lubrication  Model 

-  One  direction  of  the  fluid  flow  and  the  heat  flow  is  assumed  to 
have  negligible  effects  compared  to  other  directions 


-  The  need  for  a  mesh  in  one  direction  is  eliminated 


•  Interaction  in  the  gap  (disjoining  pressure  included)  can  be 
accounted  for  in  lubrication  model 


Using  a  Lubrication  Model  for  “inner”  region 


Particle 

Xco 

-<Z> 

Jnm 

p 

Vr  = 

cp 

YcP~ 

A 

0  1 
P 

9  Cp 

i 

cp 

V 

n 

n 

d,  u„,  Vp,  Tp 


5 


.  Vi-  T 


Transformed  NS  equations  from  x,y 


Continuity: 


dU  dV 

- 1 - 

d%  drj 


u-momentum  (v-momentum  and  Energy  equation  are  very  similar): 
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Disjoining 
Pressure 
terms  treated 
as  body  force 


Where:  u,v ->*,y  velocities; 
U,V  ->  5.  ri  velocities, 
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^Orthogonal  coordinates  (q2=0)  and  do  scaling  analysi 


^*  (Z  *71 

Scaling:  q  =—  rj  =  — 

Z  a 


By  assuming  small  curvatures  in  the  interfaces  and  using  the  fact  that: 

^  =  uyn  ~vxn  V  =  vxn  -uy # 


To  get: 

Mass: 


U-mom: 


V-mom: 


Energy: 
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Notice  that  p*  and  n*  only 
depend  on 


Pressure  Equation  in  the  gap 
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BCs:  At  5=0,  p=pouter,0 
At  5=L,  p=p0Uter,L 


Temperature  Equation: 
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Coupling  the  outer  flow  equations  to  the  inner 

lubrication  equations 


BCs:  At  £=0,  P= Pouter 

At  f=L,  p=pouter 


Edge  of  the  lubrication  layer  is  8  mesh  pts.  ->  enough  resolution  for 
outer  flow 


liquid 

Solved  via  A/ 

t  p 

Navier-Stokes 
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particle 


How  are  the  Boundary  conditions  applied  to  the 
outer  flow  (solved  via  NS)? 


BCs:  at  the  edges  of  lubrication  layer,  find  a  parabolic  velocity  profile 
that  satisfies  mass  conservation  (treated  as  outlets) 


Forces  acting  on  the  bodies 


The  forces  acting  on  the  bodies  can  now  be  obtained  from  the 
outside  flow  and  the  inner  lubrication  flow,  so: 


^  F  —  F INNER  +  FqUTER 
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The  temperature  of  the  interfaces? 


Temperature  of  the  solidification  front  is  found  via: 


T:  =  T 


m 
U  J 
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premelting 


T.  =  T. 
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No  premelting 


The  d-distribution  is  found  from  the  level  set  information 
The  particle  interface  temperature  is  found  from: 
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particle 


solid 


Flowchart  of  calculation 


Results  for  particle+stable  (non-dendritic) 

solidification  front 


Validation  of  coupled  model 


-Compare  model  solution  to  an  analytic  solution  for  a  simplified  case 
-Particle  approaching  flat  stationary  interface  (no  disjoining  pressure) 

-d/Rp  =  0.02,  Particle  velocity  =  500  microns/sec,  Rp  =  1  micron 


1.  Gary  L.  Leal,  Laminar  flow  and  convective 
transport  processes  :  scaling  principles 

and  asymptotic  analysis  (1992). 


Validation  of  coupled  model 


-Compare  model  solution  to  an  analytic  solution  for  a  simplified  case 
-Particle  approaching  flat  stationary  interface  (no  disjoining  pressure)1 

-Particle  velocity  =  500  microns/sec,  Rp  =  1  micron 

F  VS.  dmin/Rp 


1. 


0.000  0.005  0.010  0.015  0.020  0.025  0.030 

dmir/Rp 

Gary  L.  Leal,  Laminar  flow  and  convective 
transport  processes  :  scaling  principles 

and  asymptotic  analysis  (1992). 


Planar  Solidification  front  +  particle 


Solidification  velocity  =  500  microns/sec,  Rp  =  1  micron,  Hamaker  =  -8E-19  J, 

kp/k,  =  1.0  (planar),  no  premelting 


(S/UJ)  SA 


Velocity  vs.  t  and  d  vs.  t  plots  -  premelting 


Solidification  velocity  =  500  microns/sec,  Rp  =  1  micron,  Hamaker  =  -8E-19  J, 

k  /k,  =  1.0  (planar),  premelting 


k£/k,  >  1 .0  ALWAYS  ENGULFS! 


(s/ui)  SA 


What  about  kp/k,  <1.0 


Solidification  velocity  =  245  microns/sec,  Rp  =  1  micron,  Hamaker  =  -8E-19  J, 

kp  /k,  =  0.01  ,  premelting 


This  one  is  engulfed! 


(S/LU)  SA 


What  about  kp/k,  <1.0 


Solidification  velocity  =  230  microns/sec,  Rp  =  1  micron,  Hamaker  =  -8E-19  J, 

kp  /k,  =  0.01  ,  premelting 


This  one  is  Pushed! 


Vp  (m/s) 


Phase-Space  Plots  --  engulfment 


kp/k,  =  1.0  kp/k,  <  1.0 


NOTE: 

•  Similarity  to  Rempel  and  Worster  -  steady-state 

•  kp/k|  =  1 .0  Rempel  and  Woster  claim  pushing  at  certain 

values,  current  claim  is  always  engulfment 


Vp  (m/s) 


Phase-Space  Plots  -  engulfment/pushing 


kp/k,  <1.0 


kp/k,  <1.0 


min  umin 


•  2  “attractors”  in  phase-space  -  engulfment/pushing 

•  Slightest  change  causes  drastically  different  path  to  different  region 

•  Reason:  Interface  never  gets  a  chance  to  “bend” 


Conclusions  of  results  section: 


1 )  For  kp/k,  >  1 ,  engulfment  always  occurs 

2)  For  kp/k,  <  1 ,  the  critical  velocity  can  be  found  from  dynamical 
considerations  and  occurs  the  moment  the  solidification  front 
velocity  is  such  that  the  particle  “outruns”  the  tip  velocity 

Main  Contributions: 

1 .  Developed  a  multiscale  numerical  approach  to  solve  particle- 
solidification  front  problems  in  a  general  manner 

2.  Model  captures  essential  physics  that  past  research  has  neglected 

3.  First  to  look  at  the  physics  from  a  nonlinear  dynamics  point-of- 
view  (i.e.  phase-space  point  of  view) 

4.  Model  can  be  generalized  to  solve  a  variety  of  problems  involving 
interacting  interfaces 


Overall  Conclusions 


A.  Multiscale  method  was  developed 

1.  Capturing  the  interaction  dynamics  between  two  arbitrary 

solid  objects 

2.  Coupled  two  distinct  regions  (“inner”  and  “outer”)  using  a 

‘matched  numerics  approach’ 

B.  Application  of  interest 

1 .  Particle-solidification  front  interactions 

a.  Engulfment  occurs  for  k  /k,  >  1 . 

b.  For  kp/k,  <  1 ,  the  critical  velocity  is  found  using 
dynamical  considerations 

c.  No  assumptions  are  made  as  to  the  force  equations  or 
the  way  in  which  the  critical  velocity  is  obtained  (i.e.  no 
cut-off  gap  thickness  value),  etc... 


